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SUMMARY 


An experimental investigation is made of the purely 
torsional stall flutter of a two-dimensional wing pivoted about 
the midchord, and also of the bending- tors ion stall flutter of 
a two-dimensional wing pivoted about the quarterchord. For the 
purely torsional flutter case, large amplitude limit cycles 
ranging from + 11 to + 160 degrees were observed. Nondimens ional 
harmonic coefficients were extracted from the free transient 
vibration tests for amplitudes up to 80 degrees. Reasonable 
nondimensional correlation was obtained for several wing con- 
f igjarations . For the bending- tors ion flutter case, large ampli- 
tude coupled limit cycles were observed with torsional amplitudes 
as large as + 40 degrees. The torsion amplitudes first increased, 
then decreased with increasing velocity. Additionally, a small 
amplitude, predominantly torsional flutter was observed when the 
static equilibrium angle was near the stall angle. The general 
trends found here might.be of interest in related studies of 
wing stall flutter. 
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NOMENCLATURE 


A 


b 

C(k) 

^LI'^LR'^LO 
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P 

G 

g 
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I 



i 



L 

A 

M 

M 

aero 


“s 


Amplitude of oscillation 

Nondimensional location of elastic axis 

Semichord = c/2 

Theodorsen function 

Lift coefficient = L/%pV^S 

Harmonic lift coefficients 

2 

Moment coefficient = M/^spV Sc 
Harmonic moment coefficients 
Chord 

Real part of Theodorsen function 
Imaginary part of Theodorsen function 
Gravity 

Structural damping 

Bending displacement {positive up) 

Bending displacement centershift 

Moment of inertia of wing about pivot 
Virtual aerodynamic moment of inertia 

Total moment of inertia =1+1^ 

AT". 

Reduced frequency = oib/V 
Torsional spring rate 

Bending spring rate 

Torsional spring rate including gravity effects 
Lift 

Span of wing 
Moment 

Aerodynamic moment 
Spring moment 
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S 

S 

c 

t 

V 


a 


a 

a 
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X 


“o 

Aa 


?ip/ ?p 




y 

y 

p 

w 


“n 


mass of wing 

Aerodynamic moments defined by Eq. (12) 

Nondimens ional radius of gyration = ^x/mb^ 

Wing area 

Static mass unbalance of wing 
Time 

Velocity of flow 
Work due to friction 

Nondimens ional location of center of gravity = S /mb 

Angle of attack 

Centershift in angle of attack 

Angle of attack for no spring moment 

Initial angle of attack 

Static equilibrium angle 

Incremental angle 

Critical damping ratios (total and friction, 
respectively) 

Critical damping ratios for bending and torsion 
motions 

3 

Inertia density ratio = 4l/pSc 

2 

Mass density ratio = ra/7rpb i 
Air density 

Frequency of oscillation 
Natural frequency 

Natural frequencies for bending and torsion 
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METRIC CONVERSIONS 


1 cm 
1 m 

1 ra/sec 
1 kg 
1 kg-m 
1 kg~m 
1 kg/m^ 
i N 
1 N-m 
1 N/m 

2 

g = 9.82 m/sec 


.394 in 
3.28 ft 
3.28 ft/sec 
,0685 slugs 
.225 slug-ft 
.738 slug-ft^ 
.001942 slugs/ft^ 
.225 ibf 

.738 Ibf -ft 

« 

.0685 Ibf/ft 
32.2 ft/sec^ 
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SECTION 1 


INTRODUCTION ‘ 

The present investigation attempts to explore the basic 
nonlinear flutter behavior of two-dimensional wings. There is 
much information, both theoretical and experimental, on linear 
flutter behavior, i.e., the onset of small amplitude, self- 
excited flutter- Much less is known, however, about the subse- 
quent nonlinear behavior of the large amplitude flutter motions. 
For example, such a basic problem as predicting the steady-state 
amplitude of an actual wing flutter oscillation (assuming the 
wing is flexible enough so it doesn't break) is presently an 
extremely difficult task to do., , At these large amplitudes of 
motion, aerodynamic stalling occurs and alleviates the aerodynamic 
forces. Further, because of the nonlinear aerodynamic forces, 
a wing that may be stable to small disturbances, might go into 
a large amplitude flutter oscillation, if given a large enough 
disturbance. Also, structural nonlinearities may be present which 
can affect the amplitude of the motion. Such considerations of 
large amplitude flutter (stall flutter) may play important roles 
in turbomachinery, helicopters, flexible wing aircraft, and 
building structures. 

In a previous report involving one of the present authors, 
an experimental investigation was made of the purely torsional 
stall flutter and nonlinear divergence of a two-dimensional flat 
plate wing, pivoted about the midchord. See Dugundji and 
Aravaraudan^. The present report continues that work by investi- 
gating further variations of the basic wing there and attempting 
to correlate the results nondimensionally. Also, the present 
report begins an investigation into the large amplitude, bending- 
torsion stall flutter behavior of a wing section pivoted about 
the quarterchord. 

. There have been many investigations of stall flutter in 

2 . 3 

the past. The early work of Victory , Bratt and Wight , Half man, 

4 5 6 

Johnson, and Haley , Rainey ' , to mention a few names, has been 
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supplemented by more recent work on the subject (see for example. 

Refs. 7 through 15) . Typical discussions of stall flutter 

1 

applications to turbomachinery are given by Carta , to helicopters 

i 7 1 8 

are given by Crimi , and to the space shuttle is given by Reed 

It is hoped the present report will contribute some additional 

knowledge to the behavior of wings oscillating at large amplitudes. 
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SECTION 2 


TORSIONAL STALL FLUTTER OF A 
FLAT PLATE WING 


2 . 1 Experiment 

The experimental set-up for purely torsional stall flutter 
is shown in Fig. 1. A small wind tunnel with a 38.1-cm. x 5.7-cm. 
test section was used. One side wall of the test section was 
made of plexiglas for ease in viewing the test model -behavior. 

The other side wall had an easily removable wood panel, 38.1— cm. 

X 38.1 cm. square, to which the test model and its supporting 
assembly were attached. The velocity of the wind tunnel could 
be adjusted continuously from 0 to 17 m/sec. 

The two-dimensional flat plate wing was .64-cm. thick x 
10.16-cm. chord, was pivotdd. about the midchord, and had a rounded 
leading edge and a square trailing edge. The wing was restrained 
by an .074-cm. diameter steel wire acting in torsion, which pro- 
vided a linear torsional spring over very large angles. A low 
friction potentiometer measured the angular position of the wing, 
while a strain gage at the attachment block measured torsional 
moment. Also, angle lines were drawn every 5 degrees on the 
inside wall of the tunnel for visual observation of the angular 
position. The initial angle of attack could be varied through 
360 degrees by rotating the entire wing support assembly. 

Further details of the experimental set-up are given in 
Ref. 1. The overall test apparatus seemed to provide a relatively 
simple, low friction, linear torsional spring device to investigate 
large amplitude stall flutter. The nonlinearity entered here 
primarily through the aerodynamics, and no complicating structvural 
nonlinearities were present. 

In Ref. 1, the above set-up was used to take measurements 
of static moments, static divergence positions, flutter ampli- 
tudes, center shifts, frequencies, and decay-growth rates for a 
basic flat plate wing configuration. From these, nondimens ional 
harmonic coefficients, Cj^ = + i were deduced for the flat 

plate at an initial angle of attack = 0°. 
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In the present report, the runs on the basic wing configura- 
tion were extended to other initial angles = 5°, 10°, 15°, 
and the corresponding harmonic coefficients and steady 

coefficient were deduced. Additionally, runs were made on 
the wing with a large added moment of inertia which decreased 
the torsional frequency from 8.2 Hz to 4.25 Hz. Also, runs were 
made on the wing with a heavier torsional spring (.091-cm. dia- 
meter) which increased the torsional frequency from 8.2 Hz to 
12.0 Hz. Finally, the effect of nose bluntness was explored by 
pivoting the wing 180 degrees so that the square section end 
became the leading edge and the rounded end became the trailing 
edge. 

The actual procedures- for,, carrying out the above test runs 
were the same as those described in Ref. 1, and are discussed 
in more detail there. However, it should be mentioned again 
here, that, because of considerable blockage of the wind tunnel 
at high angles of attack, the velocity for any of the test runs 
was defined to be the velocity registered by the tunnel manometer 
when the wing section was brought to zero angle of attack manually 
and the variac power setting for the wind tunnel motor was not 
changed . 

Figure 2 shows the measured static aerodynamic moment 
characteristics for the flat plate. These were obtained by 
replacing the basic .074-cm. diameter torsional spring wire by 
a much thicker .318-cm. diameter wire to prevent appreciable 
elastic twist, and then using the strain gage at the attachment 
to record the torsional moment. The measurements indicate a 
linear range from’ -9 to +8 degrees and a stalled region there- 
after. 

The static divergence characteristics of the flat plate 
are shown in Figs. 3a, b, c for the basic wing case, (, 074-cm. 
torsional spring) , increased inertia case, and increased spring 
c^se (.091-cm. torsional spring). These were obtained by giving 
the wing a specified initial angle of attack a^, and then varying 
the tunnel velocity in small increments from 0 to 13 m/sec. The 
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resulting static angles were obtained visually by viewing 
the wing against the angle lines drawn on the inside wall of the 
tunnel. For some of the lower initial angles it was necessary 
to lightly hold the wing pivot from the outside to prevent self- 
excited flutter oscillations from building up around the static 
equilibrium position. The measurements show a sharp rise in the 
otQ=0 curves in the neighborhood of 5.5 m/sec for the basic wing 
and increased inertia cases, and about 8.0 m/sec for the increased 
spring case. These sharp rises represent a static divergence 
phenomenon. However, the nonlinear characteristics of the aero- 
dynamic moment keep the static equilibrium angles finite. 

The torsional stall flutter characteristics of the flat 
plate wing were investigated by giving the wing a specified angle 
of attack a^, and then varying the tunnel velocity in small incre- 
ments from 0 to 13 m/sec, as for the static divergence tests. 

At a given velocity, the wing section was manually given various 
disturbance angles Aa from the equilibrium position, then let 
go, and the resulting transient potentiometer position signal 
was recorded. This gave much data for transiently analyzing the 
response. Also, the amount of disturbance Aa needed to cause 
the wing to go into a steady flutter limit cycle oscillation was 
noted as well as the oscillation amplitude, centershift, and 
oscillation frequency. Before and after each test run, a static 
calibration was taken for the angular position signal, and also 
a transient decay record from a 90 degree twist angle for the 
natural frequency and critical damping ratio at zero velocity. 

Figure 4a shows the measured amplitudes of the stall flutter 
for the basic wing case. Some of these amplitudes reached + 100 
degrees at the higher velocities. The dotted lines indicate 
roughly the amount of disturbance Aa required to cause the stall 
flutter limit cycle. It can be seen that self-excitation occurs 
for initial angles a^ between 0 and 8 degrees. Thereafter, 
increasxngly larger disturbances Aa are required as a^ increases. 
Apparently, above a certain critical velocity, stall flutter is 
initiated if the wing angle gets near the stall angle of about 
10 degrees, either by sufficiently large disturbances Aa (for 
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high a^) or by the natural divergence process (for low a^) . 

Also, the gradual change in the flutter amplitude curves from 
a^=0° to 15° is apparent in Fig. 4a, In fact, for a^=5°, two 
different amplitude flutter limit cycles can occur at V ~ 5.5 
m/sec depending on the initial disturbance Aa, A short movie 
was made of several of these limit cycles for the basic wing 
case. 

Figures 4b and 4c shown the corresponding flutter ampli- 
tudes for the increased inertia and increased spring case respec- 
tively, The increased spring curves are similar in shape to the 
basic wing case with the velocities for a given amplitude increased 
proportionally. The increased inertia curves have a somewhat 
different shape with more double limit cycles over the range 
a^=5° to 15°, and the velocities for a given amplitude are de- 
creased from the basic wing case. 

The flutter oscillation frequency w and the centershift 
at flutter were also recorded and plotted* for the three wing 

j 

cases, but are not shown here. Generally, the flutter frequencies 
were slightly lower than the natural frequencies except at 
low flutter amplitudes A < 50°, where they decreased somewhat. 
Similarly the centershifts were close to the initial angles 
except at low amplitudes where they increased somewhat. Curves 
of flutter frequency and centershift for the basic wing case are 
shown in Ref, 1. 

As mentioned previously, many transient responses were obtained 

\ 

in addition to the steady state flutter limit cycles. The ampli- 
tudes of these transient decays (or growths) were plotted versus 
cycles (effective time) on semi-log paper, and from the resulting 
slopes of these semi-log plots, the critical damping ratios ^ 
could be obtained as functions of amplitude A. Also, the frequency 
« and center shift were obtained as functions of amplitude A. 

Such plots of 5, oj, a versus amplitude A are shown in Ref. 1 for 

o 

the basic wing case. It should be mentioned that all transient 
responses were read carefully with a magnifying glass and the 
amplitudes, frequencies, and centershifts were subsequently 
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smoothed by a numerical averaging process. 

To obtain the purely frictional damping for the basic wing, 
the increased inertia wing, and the increased spring wing cases, 
the transient decay records at zero velocity were repeated with 
the wings replaced by a small 3.18-cm. diameter steel disk having 
the same moment of inertia as the wing. These disk results pre- 
sumably had negligible aerodynamic damping and hence were indica- 
tive of the friction in the bearings and potentiometer. Values 
of the friction damping ratios Cp for the three cases are given 
in Fig. 5. The critical damping ratios increase as the amplitude 
A decreases. This is indicative of Coulomb type damping. Also, 
the for the increased spring seems somewhat lower than for 
the basic wing and increased inertia cases. 

In the analysis of the transient response records at a given 
cunnel velocity, it is well to point out that obtaining the total 
damping ratio as a function of amplitude was not always a pre- 
cise procedure, and in some cases involved some judgment in read- 
ing the response traces. This arose from the presence of "over- 
shoots" in some of the responses which had to be faired in 
reasonably. Examples of these are given in Ref. 1. Generally, 
the basic wing and increased inertia wing had fewer of these 
"overshoots" and the total damping results were reasonably 
indicative of the actual behavior. The increased spring wing, 
however, had many of these "overshoots" and the transient response 
analysis was not done for this case. Only steady-state oscillation 
cases were used here. 

2.2 Theory and Discussion 

The static aerodynamic characteristics of the flat plate wing 
about its midchord are expressed by the nondimensional coefficient 



a, 

For a two-dimensional flat plate in incompressible flow, the 
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theoretical moment curve slope in the linear (unstalled) region 
is 

_ IL =: 1,57 

doc 2_ (2) 

for the wing pivoted about the midchord. This theoretical value 
compares reasonably with the experimentally measured value of 
1.39 per radian shown in Fig. 2. The somewhat lower value here 
is mainly due to the blunt trailing edge and associated viscosity 
effects which prevent the full theoretical lift-curve slope of 
2v from being achieved. 

The linear divergence speed for the wing can be found by 
considering static moment equilibrium about the midchord, namely, 

= O (3) 

where is the torsional spring rate and is the initial angle 
of attack. Setting = (dCj^/da)a in Eq. (1), then solving for 
a, then noting when a results in the linear divergence speed, 

(4) 

The measured experimental values for the basic wing, the increased 

inertia wing, and the increased spring wing are given in Fig. 6. 

Placing these values into Eq. (4) gives calculated divergence 

speeds of 5.7 m/sec, 5.7 m/sec, and 8.3 m/sec respectively for 

the three cases. These linear divergence speeds are indicated 

in Figs. 3a, b, c by a small arrow and agree well with the sharp 

rise in the observed static angle in that vicinity for the 

a =0 cases, 
o 

The actual static divergence characteristics of the flat 
plate wing are governed by the nonlinear versus a curve shown 

in Fig. 2, rather than the linear relation, Eq. (2) , which is 
valid only for a between -9 to +8 degrees. To obtain these 
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nonlinear static divergence characteristics, one rearranges 
Eq. (3) in the form. 


c 


M 



((X- oc,) 


(5) 


The left-hand side of Eq. (5) represents the versus a curve 
of Fig. 2, while the right-hand side represents a straight line 
of slope 2k^/pV Sc which passes through the point a=a^. A 

graphical solution of Eq. (5) is easily obtained by locating the 
intersection of the straight line with the Cj^ versus a curve. 
Note that for low velocities, there will be only one solution 
while for higher velocities there will be three {but the second 
solution is always statically unstable) . 



Graphical Solution of Eq. (5) 
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Figures 3a, b, c show the calculated static equilibrium 
angles as determined by this method for the three flat plate 
wing cases. Good agreement with experiment is found over the 
given range of initial angles for all three wing cases. It 
is interesting to note that the nonlinear divergence phenomenon 
here resembles the nonlinear buckling of plates with varying 
amounts of initial imperfection. 

The dynamic behavior of the wing torsional system with no 
aerodynamic forces present can be represented by the standard 
equation, 

I + Z (X + = 0 


( 6 ) 

where I is the moment of inertia about t!\e midchord, = /k^/I 
is the natural frequency of the wing, and 5 ^, is the critical 
damping ratio due to the bearing and potentiometer friction. 
Experimental values of k^, and I for the basic wing, the 

increased inertia wing, and the increased spring wing are given 
in Fig. 6. In these and were measured directly from static 
and free vibration tests while I was determined from these 
measured values of k^ and by using the relation, I = k^/iOj^. 

The amount of virtual inertia of air that would contribute to 
this value of I, can be estimated from two-dimensional aerodynamic 
theory (Ref. 19) as. 


I 


V 


IZ8 


(7) • 
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2 

This result is I = .18 x 10 kg-m which would give a 
negligible contribution here. The friction critical damping 
ratios for the three wing cases were measured from transient 

r 

decay records as mentioned earlier, with the wing replaced by a 
solid disk of equivalent inertia. The results are given in 
Fig. 5 and are seen to be functions of oscillation amplitude A, 
thereby indicating nonlinear damping behavior. The qualitative 
form of this friction damping can be deduced roughly by consider- 
ing the concept of "equivalent damping". Ref. 20, as follows. 

For the bearing and potentiometer friction, the resisting 
moment is apt to be predominantly a Coulomb friction-type 
moment of constant magnitude, i.e., Mp=c^. Assiiming harmonic 
motion, a=A cos wt, the work done over one cycle is, 

W = 4- [ i • dt = 4 o. A 

J . (8) 

0 

For the usual viscous-type resisting moment = 2? I a, the 
work done over one cycle is, 

^ ; ■ - 
Wp = 4 1 « tit = Z t U I TTW A 

(9) 


Equating the work done for these . two cases gives the equivalent 
critical damping ratio as. 



.634 Oi 

JkrA 


( 10 ) 


since lo - (o^ here, the 5 ^, should vary inversely with amplitude A. 
Also, increasing the stiffness will lower the damping ratio Cp, 
whereas increasing the inertia will have no effect. All these 
trends are borne out by the experimental results given in Fig. 5. 
Of course, -some viscous and other type dampings may also be 
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present in the actual friction which may modify the theoretical 
results somewhat, but Eq. ( 10 ) characterizes the predominant 
components here . 

The stall flutter characteristics of the flat plate wing 
at a wind velocity V can be represented by the equation, 

Ic< +■ (o( - CKp') = M (11) 

where the aerodynamic moment is given by 

i ^ C. (ok- 0(c) Y ] “ (12) 

In the above, represents the natural frequency without the 
air, b is the semichord, is the aerodynamic virtual mass, 
while m 2 , m^^, and m^ specify the additionaH aerodynamic damping, 
spring, and constant moment, respectively, due to the wind velocity 
The aerodynamic quantities m 2 , m^, and m^ are nondimensional, and 
are assumed here to be functions of oscillation amplitude A, center- 
shift about which the oscillation takes place, and reduced 
frequency wb/V (actually, the static mp term may be considered 
roughly independent of wb/V) ^ One can then combine Eqs. (11) and 
( 12 ) into the standard form, 

2 . 

0 (+Z'ip^U) 0 ( + 600 < = CO (13) 

where 5 represents the total critical damping ratio, w the 

J. ^ 

natural frequency, and a the center shift, respectively, of the 
system response in the presence of the wind velocity. The quanti- 
ties 5 , CO, a generally vary slowly during the response, and are 
X o 

assumed here to be functions of oscillation amplitude A. 

From transient decay records of the wing at a given velocity 
V and initial angle of attack setting a^, the critical damping 

*Note, there is a slightly different definition of m^ and m^ here 
than was given in Ref. 1, because of the term. Also, vxrtual 
mass is included here. 
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ratio frequency w, and center shift a can be obtained as 

functions of amplitude A. Then from Eqs. (11) , (12 ) , (13 ) , one 

can determine the aerodynamic quantities m 2 r of Eq. (12) 

as. 


m, = 




(14) 


(15) 


(16) 


3 

where pi s 4l/pSc is an inertia density ratio parameter, while 
= I + represents the total moment of inertia of the system. 
The virtual mass can be either estimated from theory using 
Eq. (7) , or can be deduced from natural frequency measurements 
at zero velocity with and without the fluid medium present. If 
air is the fluid medium (as here) , then the tests must be done 
in air and in a vacuum; if water is the fluid medium, the tests 
can usually be done in water and in air. 

The aerodynamic moment given by Eq. (12) can also be 

expressed in the standard harmonic coefficient form by asstiming 
harmonic motion about some centershift a as, 

c cotr 

Oc(-fc) = ■ 0(<, + Ae (17) 

and then writing, 

M =4 ^ ^ 

asAo 2. j 
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Upon substitution of Eq. (17) into (12) , comparing with Eq. (18) 
and using Eqs. (14)., (15), (16), one can express the harmonic 
coefficients as 


c 





(19) 


c 


MR 




5^c*J(20) 


C 


M O 


m. 



(21) 


As in Eqs. (11) and (12), the above harmonic moment coefficients 

are similarly functions of amplitude A, centershift and 

reduced frequency cob/V, while the static term may be considered 

rpughly independent of wb/V. Note also that the nondimensional 

2 

parameter k^/pV Sc which appeared in the static divergence analysis 
Eg. (5), can be expressed alternatively as. 


[fmf 


Experimental values of k^, I and y for the basic wing, 

the increased inertia .wing, and the increased spring wing ^are 
given in Fig. 6. As mentioned previously, for these wing cases, 
the virtual mass -contribution I estimated from Eq. (7) was found 
to be negligible. Hence in all previous Eqs. (12) to (21) , one 
sets 1^=0 and I^/I=l. 

Figures 7a-d, 8a-d, 9a— e, show the experimental harmonic 
coefficients respectively for 'Z 0°, 6°, 11°, 16° 

and various values of amplitude A and reduced velocity, V/ub. 

These were obtained by applying Eqs. (19), (20), (21) with 1^=0, 
I,p/I=l to the measured transient response values of for 
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the three wing cases. Reasonable agreement for the harmonic 
coefficients was generally obtained by the three different wing 
cases, thereby bearing out the basic nondimensionalization pro- 
cess. Looking at the C ^ curves. Figs. 7a-d, it seems that the 
larger amplitude A=40°, 60°, 80° curves do not change much in 
shape or magnitude in going from centershifts a =0° to a =16°, 
but the A=10 and .20° curves do change considerably. In particular, 
the A=10° curve is always stable (negative C„^) for all V/mb at 
CL ~ 0 , but becomes unstable at a lower and lower V/wb as a 
increases. The lowest instability (positive Cj^j) generally sets 
in around V/bw ~ 2, except for ~ 16° where it becomes a little 
less at the lower amplitudes. Looking at the curves Figs. 8a- 

d, it appears again that the larger amplitude A=40°, 60°, 80° 
curves do not change much in shape and magnitude but the lower 
amplitude A=10° and 20° curves give lower magnitudes as in- 
creases. Looking at the curves Figs. 9a-e, it appears that 

despite the scatter, the coefficient correlates reasonably 

with and A only. A previous attempt to correlate additionally 
with V/buj as for and gave a meaningless jumble. 

The linear, two-dimensional theory solution for a flat plate 
oscillating in rotation about its midchord is given as (see Ref. 
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where I is the span, k = ub/V is the reduced frequency, and 

C(k) = F + iG is the Theodorsen function. Upon assuming harmonic 

motion about some center shift a as in Eq. (17) , the above gives 

c ^ 

the harmonic coefficients of Eq. (18) as. 
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This linear theory solution is also shown in Figs. 7a, 8a, 9a. 

The is always negative (always stable) and seems to agree 

reasonably with the experimental curve for A=10°. The large 

increase in C,,„ at low values of V/ojb is due to the virtual mass 
2 

{the k/8 term in Eq. (25)), and would be absent if virtual mass 
I^d were neglected. The resulting dimensioned moment remains 
small as V 0 due to multiplication of by dynamic pressure 

pV^/2. 

A theoretical estimate of the harmonic coefficients C.,_ and 

MR 

Cmo also be made based on quasi-static aerodynamic consideratior 

using the static moment characteristics shown in Fig. 2. This 

estimate would assume that under dynamic conditions, the moment 

coefficient C„ still varies with a as shown in Fig. 2, i.e., 

there is no dynamic stall overshoot. Assuming harmonic motion 

about some centershift a as in Eq. (17) , one can determine the 

c 

instantaneous as a function of time and then Fourier analyze 
M 

it to obtain the average component and the fundamental' har- 
monic component This is described in Ref. 1, and results 

are shown there for a wide variety of centershifts and ampli- 
tudes A. The appropriate quasi-static values of are shown 

in Figs. 8a-d on the right hand side. These are of course inde- 
pendent of V/o)b and hence would appear as straight lines. The 
experimental values of considerably greater than the 

quasi-static ones, thereby indicating considerably more aero- 
dynamic moment present than would be expected from simple quasi- 
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static theory. This, is probably due to dynamic stall overshoot 
effects. For large V/cob (i.e.^ slow oscillations), the quasi- 
static values appear as possible asymptotes to the experimental 
results. The appropriate quasi-static curves are shown in 
Figs. 9a-e. These experimental values of agree reasonably 

with the quasi- static ones at low amplitudes, but seem to show 
somewhat more moment at the higher amplitudes, again probably 
due to dynamic overshoot effects. 

The general inoment coefficient curves given 

in Figs- 7a-d, 8a-d, 9a-e define the basic nondimensional charac 

teristics of stall aerodynamics for this flat plate wing. From 

these C„T-, C.,_, curves, one can then obtain either the tran- 

MI MR MO 

sient decay-growth characteristics or the steady limit cycle 
flutter oscillations for other physical situations — different 
sizes, inertias, springs, and dampings. To obtain the transient 
decay - growth characteristics, one rearranges Eqs. (19), (20), 
(21) to the form. 
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One first uses Eq. (7) to estimate the virtual -Inertia I^. If 
I^/I is small (less than say 5%) one may set I^— 0 in the above. 
For a given amplitude A, centershift a , and reduced velocity 


24 



V/wb, one can determine C„„, C„^. Then from Eqs. (27)- 

(30) one finds tOjj/to, and V. By varying V/cob, one can 

obtain a plot of versus V for a given amplitude A. To obtain 
the stall flutter limit cycle, one continues to vary V/wb in 
the above process until one obtains This then gives the 

velocity V for the given limit cycle amplitude A. It can be 
seen from Eq. (28) that for no friction damping present, 
the flutter point will occur at that V/wb where first becomes 

positive . 

The above procedure was actually carried out using the basic 
^MI ' ^MR' ^MO reproduce reasonably the basic wing, the 

increased inertia, and the increased spring curves given in Figs. 
4a, b, c. In doing this, use was made of the basic wing data 
given in Figs. 5 and 6 . Additionally, calculations were made for 
the three wing cases with no friction damping present, The 

r 

results are shown in Pig. 10. For no friction damping ^^=0, in- 
creasing the spring rate shifts the velocities proportionally to 
the natural frequency u)j^. Increasing the inertia also shifts the 
velocities roughly proportionally to oj^^, but there may be some 
slight shape change due to the inertia ratio y which affects the 

value. The addition of some friction can change the re- 
sults considerably. The increased inertia case seems to be most 
affected by the addition of friction 5 ^. 

In concluding this section, it should be mentioned that the 

basic curves can be obtained either from free tran- 

MI MR MO 

sient response tests as here, or from measured force or pr^essure 
measurement on an externally forced oscillating wing as in Refs. 
11 and 6 . The forced method allows more flexibility in covering 
the required range of parameters. However, it might be argued 
philosophically that the coefficients deduced from the free 
transient response tests,- particularly when as in a steady 

flutter oscillation, are mo^e physically direct for predicting 
flutter, since they are deduced from an actually observed flutter 
situation. 


25 



2.3 Effect of Square Leading Edge 

To explore the effect of nose bluntness on the stall flutter 
and static divergence behavior, the flats plate section was 
pivoted 180 degrees so that the square section faced the stream 
and the rounded section became the trailing edge. 

Test results for the static moment characteristics are 
shown in Pig. 11 along with the previous basic wing for compari- 
son. The maximum C.. achieved now was somewhat less than for the 

M 

basic wing and the linear moment curve slope was 1.07 per radian 
instead of the previous value of 1.39 per radian. These were 
probably due to the earlier stalling due to the sharper leading 
edge and the poorer realization of the Kutta condition due to 
the rounded trailing edge. The corresponding static divergence 
characteristics are shown in Fig. 12. The experimental values 
here are roughly similar to that of the basic wing case Fig. 3a, 
although the points seem slightly lower. The agreement with 

theory is also somewhat poorer here, indicating that either some 
of the measured angles were a little off, or perhaps the measured 
curve of Fig. 11 was a little low. 

Figure 1 3a shows the measured steady-state amplitudes of 
stall flutter. Comparing these with the previous basic case. 

Fig. 4a, reveals the same general characteristics for the two 
cases. The lower amplitudes of the a^=5° case do seem to occur 
at a slightly lower velocity than the basic wing case, but the 
same two different amplitude limit cycles occur here as before. 
Figure 13b shows the flutter amplitudes for larger initial angles 
of attack, a^=20^ to 90^, for general interest. Beyond 
no limit cycle amplitudes smaller than 90° could be found, 
although larger ones were present. For a^~90°, two different 
amplitude limit cycle branches occur, one with amplitudes around 
100° and another around 145°. To obtain these, large disturbances 
Aa > 90° had to be’ given to the wing. incidentally, similar 
behavior as above occurs for the basic wing, although this was 
not reported fully in Ref. 1. 
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In suininary then, it appears that the effect of nose blunt- 
ness on the stall flutter and static divergence behavior .is 
relatively minor for this flat plate. This- is probably to be 
expected for the relatively low thickness ration, 6.3%, occurr- 
ing here. For thicker flat plates, the effect may be greater. 
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SECTION 3 


BENDING-TORSION STALL FLUTTER 
OF A WING SECTION 


3 . 1 Experiment 

A preliminary investigation was made into the bending- 
torsion flutter of a wing section. The experimental set-up is 
shown in Pig. 14. The same small wind tunnel was used here as ' 
for the purely torsional stall flutter investigation of Section 2. 
The easily removable wood panel, '38.1-cm. x 38.1 cm. square, was 
replaced by another such wood panel which now had the bending- 
torsion supporting assembly attached to it. This supporting 
assembly consisted of a vertical rod which slid up and down 
between two bearings 10,2 cm. apart to provide the translational 
motion, and a rotation bearing in the center of the rod to pro- 
vide the rotational motion. The translation motion was restrained 
by two linear springs, and the rotational motion was restrained 
by an ,066-cm. diameter steel wire acting in torsion. Both the 
translational motion and the rotational motion could be locked 
out independently by an appropriate set screw. The vertical 
position was measured by a strain gage at the base of a small, 
very flexible cantilever steel beam whose free end was attached 
near the upper end of one of the linear springs (where the linear 
spring motion was small) . The angular position was measured by 
a potentiometer. Also, vertical lines as well as angular lines 
were drawn on the inside wall of the tunnel for visual observation 
of the vertical and angular position. The initial angle of attack 
could easily be varied by changing the clamping position of the 
.066-cm. diameter restraining torsional wire. A picture of this 
supporting assembly is shown in Fig. 14b. 

The two-dimensional wing section was a roughly symmetrical 
0012 airfoil section with a 12.70-cm. chord and a 1.52-cm. thick- 
ness, which was pivoted about the guarterchord . A small weight 
was placed in the wing section at about the midchord to give a 
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substantial static unbalance. The wing had a small rod at its 
quarterchord by which it could be attached to the supporting 
assembly. In doing this attaching, a slot had to be provided 
in the tunnel wall to allow the wing to translate vertically. 

To minimize the leakage around the slot, a rubber membrane with 
a thin slit was placed over the slot to close the gap around 
the moving rod. Also, a thin aluminum disk about 13.0-cm. in 
diameter was placed at the base of the airfoil to help prevent 
pressure leaking from one side of the airfoil to the other through 
the tunnel slot. Measured values of the wing section and its 
static and dynamic properties are given in Fig. 15. 

Static tests were first carried out to measure the aero- 
dynamic lift and moment acting on the wing section. The .066-cm. 
diameter torsional spring wire was replaced by a much thicker 
.318-cm. diameter wire and an appropriate strain gage to measure 
torsional moment, while the vertical position strain gage on the 
restraining linear spring was used to measure the spring deflec- 
tion and hence lift force on the wing. The angle of attack was 
varied in small steps and readings were always taken with wind 
on and wind off to measure the resulting increments in force 
and moment due to the aerodynamics. As mentioned in Section 2, 
because of blockage in the tunnel, the velocity for this and any 
of the subsequent test runs was always defined as the velocity 
registered by the tunnel manometer when the wing section was 
brought to zero manually and the variac power setting for the 
wind tunnel motor was not changed. 

Figures 16 and 17 show the measured static lift and moment 
characteristics of the airfoil section. The lift curve shows a 
linear lift curve slope of 5.35 per radian, which agrees reason- 
ably with the theoretical value of 6.28 per radian, and stalling 
occurring at about +14 degrees and -13 degrees. The moment 
curve shows an essentially zero moment in the same unstalled 
r.egion and a large negative restoring moment outside it. There 
is a substantial discontinuity- occurring at the stalling points. 
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The zero moment in the linear range again agrees with theory 
for a symmetrical wing section pivoted about the quarter chord. 

It should be mentioned that before the pressure leakage around 
the slot was minimized by the thin alxaminum disk referred to 
earlier, the lift curve measurements indicated somewhat of a 
discontinuity in angle around a=0 degrees and a much lower 
average lift curve slope of about 2.5 per radian. 

The spring rates for the translation and rotation restrain- 
ing springs were measured next. The translation spring was 
found to be linear over its entire range with a spring rate of 
139 N/m, while the torsional spring had the nonlinear character- 
istics shown in Fig. 18. These nonlinear characteristics are 
due to gravity effects of the static unbalance of the wing sec- 
tion and will be discussed in the next section. When the wing 
was oriented vertically instead of horizontally so that no gravity 
moment acted, the measured torsional spring was linear over its 
entire range with a spring rate of .0137 N-m/rad. 

The static equilibium characteristics of the wing section 
are shown in Fig. 19. These were obtained by first locking out 
the translational motion of wing, then giving the wing a specified 
initial angle of attack a^, and then observing visually the re- 
sulting static angles a as the tunnel velocity was increased in 

o 

small increments. The measurements generally show a decrease in 
angle of attack as the velocity is increased, which reflects the 
stable moment curve characteristics given by Fig. 17. It should 
be mentioned that when the static angles came in the neighbor- 
hood of about 14 degrees (the stall region) , small amplitude 
torsional oscillations developed. These will be discussed again 
later. 

The bending-torsion stall flutter characteristics of the 
wing section were investigated by giving the wing a specified 
angle of attack a^, and then varying the tunnel velocity in 
small .increments from 0 to 15 m/sec. Before each run, static 
calibrations of the translational and angular positions were 
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taken, and also dynamic transient decay responses were taken of 
the translational motion with the angular motion locked out, and 
of the angular motion with translation locked out. From, these 
dynamic tests at zero wind velocity, the natural frequencies in 
bending w, and torsion w were obtained, as well as the critical 
damping ratios c, and z for the wing section. These dampings 

il & 

here represent bearing friction plus air damping and were rela- 
tively high here, about ~ *15 and ~ "13 respectively. 

Because of the preliminary nature of this investigation, a care- 
ful attempt at extracting the purely bearing friction was. not 
made here, as was done for the purely torsional flutter of Sec- 
tion 2. In subsequent investigations such attempts would be made 
along with a general reduction of the high friction levels found 
in these tests. 

At each given wind velocity during a test run, the wing 
section was given a smhll disturbance Ah in the translation posi- 
tion and then let go to see if flutter would ensue. In general, 
the flutter was not seif-excited because of the large friction 
present, except at the higher velocities. Figure 20 shows a 
typical photograph of the wing section during flutter. The 
large bending and torsion amplitudes obtained here are well into 
the stalling range of the airfoil. Figures 21a and 21b show the 
measured bending and torsion amplitudes for initial angles of 
attack ot^=0° and 20° respectively. Aside from the large oscilla- 
tion amplitudes here, it is interesting to notice that the torsion 
amplitude decreases with velocity after having achieved a 'peak 
value. This changes the general appearance of the flutter oscilla- 
tion with velocity. Also, it appears here that the flutter started 
in at a slightly higher velocity at a^=20° than it did at a^=0°. 

In addition to this large amplitude bending- tors ion stall 
flutter, a small amplitude predominantly torsional flutter was 
observed when the wing section was rotated such that the resulting 
static, equilibrium angle a was about 14 degrees. This "small 
stall" flutter had an amplitude of about 3.5 degrees and existed 
even if the bending degree of freedom was locked out. It was 
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evidently related to the stalling of the wing section and the 
associated sharp change in See Fig. 17. When bending motion 

was permitted as well, the torsional amplitude became more regular 
and somewhat larger. The flutter amplitudes for this "small stall" 
flutter are shown in Fig. 21c, with and without the bending degree 
of freedom. 

The flutter frequencies for the large amplitude bending- 
torsion flutter and the "small stall" flutter are shown in Fig. 22. 
The large amplitude flutter frequencies vary from 4.0 to 4.7 H 

z 

which is somewhat above the natural frequencies of to, = 3.81 and 

h 

= 2.55 Hz. The "small stall" flutter frequencies varied from 
4.5 to 7.8 Hz' which was considerably above the torsional natural 
frequency = 2.55 Hz. This was due to the stiffening effect of 
the aerodynamic moment at the stall region, as will be shown later. 

For interest, the "small stall" flutter was ihvestiaged 
further by varying the restraining torsional spring of the wing 
section. Figure 23 shows the resulting flutter frequencies. As 
expected, increasing o)^ increased the resulting flutter frequencies. 

3.2 Theory and Discussion 

The static aerodynamic characteristics of the airfoil wing 
section about its quarterchord are expressed by the nondimensional 
lift and moment coefficients as 


'■ " -k?V"3' 

M 


(.31) 


(32) 


For a two-dimensional symmetrical wing section pivoted about the 
quarterchord, the theoretical lift and moment curve slopes in 
the linear (unstalled) region are dC^^/da = 6.28 per radian and 
dCj^/da. = 0 respectively. These agree reasonably with the measured 
values given in Figs. 16 and 17. 
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The static equilibriiim positions in torsion can be evaluated 
from the nonlinear moment equilibrixim equation, 

2. 




(33) 


where Cj^ is the nonlinear aerodynamic moment coefficient shown 

in Fig. 17 and M is the nonlinear spring moment shown in Fig. 

s 

18. Equation’ (33) can be solved by trial and error for the re- 
sulting equilibrium angle a . For the present investigation, 
the experimental results for given in Fig. 19 were used to 
reconstruct the nonlinear versus a curve given in Fig. 17, by 
simply solving Eq. (33) for Good agreement with the measured 

curve was obtained, thereby providing a cross-check on the 
measured values. 

The stall flutter characteristics of the wing section oscillat- 
ing through large angles, can be represented by the equations of 
motion, 

YY) Vt - S (X ce<L 0( 

~ h + I c< 

where m and I represent the total mass and moment of inertia, 
respectively, of the wing section plus associated moving support 
assembly parts, S^,is the mass static unbalance of the wing, k^^ 
is the bending spring rate, k^ is the torsion spring rate,'’ and 
is the angular position for no torsion spring moment. The 
above equations were obtained by applying Lagrange’s equations 
to the wing section. The gravity term in Eq. (34) simply repre- 
sents a constant displacement shift and can be discarded. How- 
ever, the gravity term in Eq. (35) should be included in the left 
hand side, and hence, Eq. (35) can be rewritten as. 
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where 

as. 


represents a nonlinear spring restoring force defined 



(37) 


For values of gS^/k^ close to 1, a significant nonlinearity 
develops in the spring rate. In fact, for gS^/k^ > .1 there will 
be a buckling snap- through for the torsional spring- For the 
present case, gS^/k^ = .710 and only a mild nonlinearity develops 
The measured nonlinear spring moment M versus a is shown in Fig. 

o 

18 and it agreed well with the theoretical prediction of Eq. (37) 
If one assumes harmonic motion for the wing section of the 


form. 
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then the aerodynamic lift and moment appearing in Eqs. (34) and 
(36) can be expressed in harmonic form as. 
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where the six 

harmonic coefficients C_„, C , etc, are 

ijU J_jK 

assumed to 


be functions of amplitudes .A, , A , center shift a , and reduced 

Xa CC w 

frequency o)b/V. As an approximation, the six coefficients can be 
further divided into twelve by assuming independence of the h and 
a effects, namely, 

eXc. 
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where the assumed to depend only on Aj^, 

“c' wb/V while the .... are assumed to depend 

only on A , a and ojb/V. 
oc c 

Placing Egs. (38) and (39) into the flutter Eqs. (34) and 
(36) and using a harmonic balance scheme, one can obtain a set 
of six nonlinear equations which can be solved by some Newton- 
Raphson type iteration scheme to obtain the amplitudes and 
frequencies of the nonlinear stall flutter. This, however, was 
not pursued further in this preliminary investigation. 

A standard, linear flutter analysis of this typical section 
was performed to obtain an idea of its basic, small amplitude 
flutter characteristics. Following’ Ref. 19, and using the 
measured properties indicated in Fig. 15, the following standard 
nondimensional flutter parameters were computed as,* 
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These were then placed into the flutter determinant (for a = -.5), 
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The high value of jl reflects the associated assembly support 
equipment moving with the wing. 
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where L, , L , M are standard aerodynamic coefficients which 
h a a 

are functions of reduced frequency k = ub/V/ and the Z = (1+ig) 

2 

(oi /w) is the unknown complex eigenvalue. A standard V-g 
ot 

method of solution was used to solve for the flutter speed and 
frequency. These calculations resulted in the V-g diagram shown 
in Fig. 24, which gives a zero damping g=0 flutter speed of 
Vp = 5.2 m/sec and a flutter frequency ojp = 3.9 Hz. However, 
because of the slowly rising V-g curve for this configuration, 
the effect of structural damping can be pronounced. For the 
actual wing here, the measured critical damping ratio at zero 
velocity was found to be about ~ ~ given in Fig. 15. 

This includes both the bearing function and the aerodynamic damp- 
ing at zero velocity which was not separated as in the previous 
torsional investigation of Section 2. Assuming that bearing 
friction contributes a large part of this damping, one can 
arbitrarily estimate the friction damping as Cp ~ -10. Then, 
using the rough relationship, 

^ ^ (43) 

one obtains a rough estimate of about g ~ .20 for the structural 
damping here. From Pig. 24, this would indicate a flutter speed 
of V_ = 10.6 m/sec and a flutter frequency of o)_ = 3.4 Hz, which 

r "C 

agree reasonably with the observed flutter results. Presumably, 
if the friction level were lower, the flutter would have occurred 
at a lower velocity, and would have been self-excited. 

In contrast to the large amplitude coupled bending-torsion 
flutter, the "small stall" flutter mentioned earlier seems to 
involve only the torsion motion of the wing. This seems to be 
associated with the sharp change in moment at the stalling point 
and the associated hysterysis in the development and drop-off 
of the moment. See Fig. 17. See also Ref. 15 for some further 
details of this phenomenon. The increase in torsional frequency 
far above the zero velocity natural frequency co^ can be related 
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to the stiffening effect of the jump in moment, by considering 


small oscillations of 

the 

wing 

about the stall 

position. 
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The total moment 

here 

can 

be idealized roughly 

as. 
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where AM represents the jump ih aerodynamic moment at stall. 
The above moment M, during an oscillation, can be replaced by 
its fundamental harmonic component as, 

[VI 'S; C| oJ’t (46) 


where the Fourier coefficient can be obtained by multiplying 
Eq. (46) by sin wt, integrating over a half cycle, and then sub- 
stituting Eqs.’ (45) and (44) to give, 
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dt 


^ AM + 

Tf ^ 


(47) 


The torsional vibrations of the wing are governed by the equation 
of motion, 

I w + M = o 

(48) 

Substituting Eqs. (44) and (46) into the above, solving for the 

2 

frequency o) and noting that to = gives finally. 
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CJsing the basic parameters of Fig. 15 and taking AC ~ .09 and 

M 

A - 3.5 degrees, the above Eq. (49) gives a reasonable estimate 
of the "small stall" flutter frequency, as shown in Fig. 22. 
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SECTION 4 


CONCLUSIONS 

The investigation of the torsional stall flutter and nonlinear 
divergence behavior of a flat plate pivoted about its midchord, 
revealed the following observations. 

1. The nonlinear static divergence equilibrium properties 
were accurately predicted from the measured static 
moment characteristics of the plate. Good nondimensional 
correlation was obtained for several -wing configurations. 
Above the linear divergence speed, the .static divergence 
characteristics resembled the buckling of a flat plate 
with initial imperfections. 

2. Large amplitude limit cycles ranging from + 11 to + 160 
degrees were observed. In some cases, two different 
limit cycles were found for the same velocity. 

3. Stall flutter occurred above a critical value of around 
V/boo ~ 2. It seemed to be initiated if the wing angle 
of attack came near the stall angle, either by a suffi- 
ciently large disturbance (for high initial angles a^} 
or by the natural divergence process (for low a^) . 

4. Nondimensional harmonic coefficients C,,_, C,._ were 

MI MR MO 

extracted from the free, transient vibration tests for 
amplitudes up to 80 degrees. Reasonable nondimensional 
correlation was obtained for several wing configurations. 
From these basic curves, one could obtain the stqady- 
state and transient behavior of other spring, inertia, 
and damping configurations. 

5. The effect of nose bluntness on the static divergence 
and stall flutter behavior was found to be relatively 
minor for the 6.3% thickness ratio flate plate investi- 
gated here. 

The preliminary investigation of the bending- tors ion stall 
flutter and nonlinear , static divergence behavior of an airfoil 
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wing section pivoted about its quarter chord, revealed the 
following observations. 

1. The nonlinear static equilibrium angles could be well 
predicted from the measured static moment characteristics 
of the wing. 

2. Care had to be taken to prevent pressure leakage from 
one side of the wing to the other through the wing hold- 
ing mechanism. 

3. Large amplitude coupled limit cycles, were observed, with 
torsion amplitudes as large as + 40 degrees. The torsion 
amplitudes first increased, then decreased with increas- 
ing velocity. The occurrence of coupled mode flutter 
was reasonably predicted by linear theory. 

4. Because of large friction present, the stall flutter was 
not self -excited, except at the higher velocities. Gener- 
ally, a sufficiently large disturbance was necessary to 
initiate it. 

5‘. In addition to the large amplitude bending-torsion 
flutter, a small .amplitude, predominantly torsional 
flutter was observed when the static equilibrium angle 
of the wing was at the stall point. The amplitude here 
was about + 3.5 degrees , and the frequency was much higher 
than the natural frequency at zero speed due to the stiff- 
ening effect of the aerodynamic moment at the stall. 

The above preliminary investigation of bending- torsion could 
be extended further to explore the phenomenon in more depth. In 
future work, more care should be taken to lower the high friction 
levels present, and to lower the associated assembly support mov- 
ing mass in order to make the density ratios y more typical of 
actual construction. Also, the bending to torsion frequency 
ratio should probably be lower. 

Because of tunnel blockage, the absolute velocities and 
coefficients for the torsional flutter and the bending-torsion 
flutter tests here may not be accurate, but the general trends and 
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phenomena should be similar on other wing sections. It is hoped 
the present report has contributed some further information and 
understanding to the interesting and complicated phenomenon of 
nonlinear stall flutter. 
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Experimental Set-up for Torsional Stall Flutter 








Fig, 2 Static FIoment for Flat Plate Wing 
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16. 3 c Static Divergence ( Increased Spring ) 
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Fig. 4b Flutter Amplitudes ( Increased Inertia ) 
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Fig, 4c Flutter Amplitudes ( Increased Spring ) 
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Measured Values of Frjction Damping Ratio 
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Parameters for Flat Plate Wing ( Torsional Flutter ) 
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Fig, 9a Experimental Curves ( A = 10'" ) 
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Fig, 9c Experimental Cj^q Curves ( A = ) 
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Fig» 10 Calculated Flutter Amplitudes for No Friction 
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Fig, 11 Static Moment ( Square L.E. Wing ) 
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Fig. 13a Flutter Amplitudes ( Square L.E. Wing ) 
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Fig, 18 ' Spring Moment vs, Wing Angle of Attack 
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Fig. 19 Static Equilibrium for Bending-Torsion Wing Section 
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Fig. 22 Flutter Frequencies for Bending-Torsion Wing 
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Fig. 23 Frequencies for Small Stall Flutter ( Pure Torsion ) 
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